home *** CD-ROM | disk | FTP | other *** search
- 10 OPEN "I",#1,"CDATA" :REM Compressed data.
- 20 INPUT#1,N,L :REM Get count, sampling inverval.
- 30 K=INT(N/L):DIM B(K)
- 40 FOR I=0 TO K :REM Get compressed data.
- 50 INPUT#1,B(I)
- 60 NEXT I
- 70 CLOSE
- 80 OPEN "O",#1,"RDATA" :REM Create reconstructed data.
- 90 PRINT#1,N :REM Write data count.
- 100 W=3.141593/L
- 110 FOR I=0 TO N :REM Reconstruction
- 120 IF I MOD L = 0 GOTO 190 :REM Branch at sampled values.
- 130 G=0
- 140 FOR J=0 TO K :REM The Nyquist sum.
- 150 M=W*(I-J*L)
- 160 G=G+B(J)*SIN(M)/M
- 170 NEXT J
- 180 GOTO 200 :REM Sum done; store this value.
- 190 G=B(I/L)
- 200 PRINT#1,G :REM Write reconstructed value to file.
- 210 NEXT I :REM Go reconstruct next value.
- 220 CLOSE :REM Done
- 230 PRINT "The reconstructed data file is RDATA"
-